High frequency polarization switching of a thin ferroelectric film 
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We consider both experimentally and analytically the transient oscillatory process that arises when 
a rapid change in voltage is applied to a Baj:Sr\-xTi02, ferroelectric thin film deposited on an MgQ 
substrate. High frequency (« 10*rad/s) polarization oscillations are observed in the ferroelectric 
sample. These can be understood using a simple field-polarization model. In particular we obtain 
analytic expressions for the oscillation frequency and the decay time of the polarization fluctuation 
in terms of the material parameters. These estimations agree well with the experimental results. 
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I. INTRODUCTION 



Unique intrinsic properties make ferroelectric materi- 
als attractive both for fundamental research and ap- 
plications in devices using electro-optical, piezoelec- 
tric and other effects. Bulk ferroelectric materials, 
particularly those based on barium-strontium titanate 
BoxSri^xTiOz (BST) compounds ^, are attractive 
for high-power applications because of their high dielec- 
tric permittivity and small losses. Ferroelectric-based 
devices include ultra-fast electrically-controlled phase 
shifters for amplitude and phase control. It has been 
shown that the dielectric permittivity of a ferroelectric 
can be altered by applying an electric field. Therefore 
a fast ferroelectric phase shifter controlled by an electric 
field bias is being investigated [3] to be used for appli- 
cations in particle accelerators. Ferroelectric materials 
provide significant benefits for several applications such 
as switching and control elements. These are able to 
handle high peak and average power while maintaining 
a very short response time of less thanjust few nanosec- 
onds. According to some estimations [if, 0] the response 
time to an applied external electric field is about 10~^^ s 
for crystalline and 10~^^ s for ceramic compounds. The 
high permittivity (over 1000) of ferroelectrics makes them 
potential candidates to replace silicon oxide dielectrics 
as storage capacitors for memory devices. A solid so- 
lution barium-strontium titanate BaxSri^xTiO^ (BST) 
has a high permittivity and a composition dependent 
Curie temperature T^ which varies in a range of 30 - 400 
K. This strong dependence of the dielectric constant on 
electric field offers the opportunity to use ferroelectrics 



as tunable devices. Ferroelectric based phase shifters or- 
ganised in an array have the advantage of being cheap 
and consuming a reduced power while continuously tun- 
ing the phase of a high power microwave signal. This 
rapid electrical steering is realized by adjustirig the bias 
voltages on each element. See for example {5| and 
where the phase shifting elements based on BST thin 
film capacitors are discussed. There is a certain advan- 
tage in using elements build on thin ferroelectric films 
because of their compactness and parameter tunability. 
The presence of internal stresses in the thin film ferro- 
electrics changes their electromechanical and dielectric 
properties drastically. Because of this the dielectric per- 
mittivity in a thin film is reduced by about an order of 
magnitude compared to one for the bulk media. However 
choosing the appropriate substrate allows to adjust the 
internal stress level and tune the physical properties of 
the thin ferroelectric films. Thin films of BST deposited 
by the sputtering technique are discussed in@, 0]in ref- 
erence to the production of compact tunable capacitors. 
These elements are attractive for applications in adaptive 
impedance matching networks and tunable filters. Size 
effects in thin ferroelectric films are discussed in It 
is reported that no critical film thickness is required to 
obtain polarization switching. 

A noteworthy transient effect occurs in ferroelectric 
material subjected to an alternating electric field which 
causes polarization switching between two stationary 
states of the ferroelectric. The switching process is fol- 
lowed by high-frequency polarization oscillations around 
their stationary states. Basically the polarization dynam- 
ics of a ferroelectric near its steady state can be viewed 
in terms of a damped oscillator with an eigenfrequency 
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determined by the material parameters. An alternating 
electric field serves as the external force that pushes the 
oscillator away from its equilibrium. The generation of 
infrared (IR) radiation by means of polarization switch- 
ing was proposed for ferroelectrics Q . There the authors 
estimate the energy radiated using the dipole approxima- 
tion. 

The study of the transient dynamics of these processes 
can help understand how the system parameters can be 
adjusted to provide the required transient behavior. That 
is important for devices subjected to sharp/shock peri- 
odic or aperiodic forces of high frequency. They then 
spend most of their time in a transient state, even if 
the relaxation time is smaller than the observation pe- 
riod. Transient processes occurring in ferroelectric might 
become unwanted effects if the possible applications re- 
quire fast switching between polarization states, such as 
in FeRAM. In opposite they might be required if the po- 
larization switching is used to produce oscillations that 
will form an IR pulse. The polarization relaxation time 
and oscillation frequency are important characteristics of 
this transient behavior. The polarization damping con- 
stant in ferroelectrics is of the order of lO^'^s"^, as indi- 
cated in [l3j . The theoretical description for an isotropic 
paraelectric in the framework of the dynamical Landau- 
Khalatnikov model being investigated in [ll|, [l^ allows 
to calculate nonlinear susceptibility coefficients. 

Here we discuss an experiment where we observe the 
polarization dynamics in a thin EST film. A similar ex- 
periment studying fast polarization switching was per- 
formed earlier by some of the authors. The method of 
observation of the polarization uses the second harmonic 
generation in the thin film of the EST solid solution. It 
is discussed in In the present experiment we follow 
a pump-probe procedure. First we apply a constant field 
Es and obtain the steady state induced polarization Pg. 
Then we send in an additional small electric field pulse 
to probe the film. In principle this time dependent solu- 
tion could be described by the theory set up in [l^ but 
we chose in this first study to describe the relaxation re- 
sponse of the film. For that we set up a scattering theory 
formalism. The bound states of the system {E, P) will 
then give the response of the film, ie the typical oscilla- 
tion frequency and the radiative decay time. We identify 
two channels of damping, the radiative damping and the 
inner damping. From the experimental data we estimate 
the magnitude of both terms and find that here the inner 
damping dominates. There could be other experiments 
where the contrary happens ie the radiative damping may 
dominate. 



II. MODEL OF ELECTROMAGNETIC 
RESPONSE IN FERROELECTRIC FILM 

The theory of ferroelectric has started to develop in the 
1930s. T he p hase transition theory proposed by Lan- 
dau [iHl, jl6l | was applied to describe the behavior of 



the ferroelectric near a critical point of the phase transi- 
tion. Following the Landau-Ginzburg-Devonshire theory 
the thermodynamic Gibbs potential $ in the neighbor- 
hood of the critical point can be represented as a power 
series of the order parameter, namely the polarization 
[13, m, O up to 4-th order, and ^20] up to terms of 
order P^. 

$ = $0 + + - EP. 

The expansion coefficients are a = ao{T — Tc) and 
/3. Tc is the Curie temperature. In a solid solution 
BaxSri^xTiOz the Curie temperature depends on the 
relative concentration x of barium. The approximate 
formula to calculate in EST was proposed by [2l| . 
It reads = SGOx + 40. Particularly = 2mK for 
X = 0.7. This is the concentration of barium in the thin 
film investigated in an experiment on ferroelectric switch- 
ing[l|. 

We describe an experiment where a thin film of thick- 
ness I is deposited on a dispersionless substrate. The 
optical properties of the surrounding media can be char- 
acterized by the refractive index n[z), such that n[z) = 1 
for z < and n{z) = n > 1 for z > 0, i.e. in the 
substrate. The electric field E{z, t) is governed by the 
Maxwell equation and the polarization obeys the equa- 
tion of a damped oscillator driven by the field (compare 
to [13]) 

d'^E n^d'^E I d^P,, , 

-d{z), 

(1) 



9z2 



d^P 



dP 
eo dt 



aP + pP^ ^Ef{t), 



eo dt^ 

where Ejit) is the electric field inside the thin film. The 
constant r is proportional to the inverse of the Eorn fre- 
quency, this value will be estimated for EST later. Typ- 
ically the decay constant 1/(7t2) is small [l^. However 
we take it into account for generality. 

The electromagnetic wave is incident from the left of a 
film located at z = 0. Thus the electric fields outside the 
thin film E~{z,t) for z < and E'^{z,t) for z > are 
defined by free wave equations. The time Fourier images 
of these fields can be written as 



E-{z,u) = Ae''''' + Be-'''''', 

E+{z,uj) = Ce'''^'. 

Here ki = uj/c and = urn/c are the wave numbers on 
the left and on the right of the thin film, A is the Fourier 
amplitude of the incident wave, B is the Fourier ampli- 
tude of the reflected wave and C is the Fourier amplitude 
of the transmitted wave. 

According to the boundary conditions at z = [l^, 
[2^ the electric field E and its spatial derivative E,z are 
connected by the relations 

M-^z = 0,uj) = E^+\z = 0,u}), 

h ? 

E\-\z = 0,u;) - £;(+)(z = 0,w) = 4^P(a;), 
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FIG. 1: Schematics of the experimental setup where an elec- 
tric field pulse polarizes a ferroelectric thin film deposited on 
a substrate. 



where P{uj) is the Fourier image of the thin film polariza- 
tion. This leads to the following equations for the Fourier 
amplitudes of the left and right electric fields 



A + B ^C, A-B = nC 



iluj 



ceo 



Pico) 



Hence the amplitude of the transmitted wave C and 
the amplitude of the refiected wave B can be expressed 
through the amplitude A of the incident wave and the 
thin film polarization P{uj): 



C = 



-A 



ilu! 



n+1 c{n + l)eo 



B 



1 -n 



-A- 



ilu) 



l + n c{n + l)eo 



P{i0) 



(2) 
(3) 



Using the inverse Fourier transform we obtain the am- 
plitude of the electric field inside the thin film Ef as 
E+{z,t) at z = [Hj, i.e., 



Efit) 



l + n 



I 



dP 



c{n + l)eo dt 



where Ein{t) is the electric field of the incident wave. 
Thus one can find from the second equation of the system 
((!]) that the polarization of the ferroelectric thin film is 
governed by the following equation 



^8FP_ 

eo dH 



7T 



aP + PP'' 



dP 

c{n + l)eo J dt 
2 



l + n 



E^nit). 



(4) 



This expression shows that two relaxation channels ex- 
ist. One of them is the ordinary one related with inner 
friction i.e., the 7-term. The second relaxation channel 
is due to the radiation process. The variation of the po- 
larization generates the electromagnetic field outside the 
film. 



The polarization of the thin film evolves according to 
(HI as an oscillator with damping and a forcing that is 
equal to Eg = 2Ein{t)/{l + n). Furthermore, this equa- 
tion allows to consider as a constant electric field i.e. 
a constant voltage. In this case Eg — Vg/l where Vg is 
the applied voltage. 



A. Relaxation to a steady state polarization 

In the absence of an external electric field the fer- 
roelectric possesses two equilibrium states, each corre- 
sponding to different polarities. In the paraelectric phase 
there is only one equilibrium corresponding to an unpo- 
larized state. The experiment under consideration was 
performed at room temperature which is above the Curie 
temperature for EST so that the sample is in the para- 
electric state. However there is another way to change 
the polarization state. This can be induced by applying 
an external field or a stress. In this experiment we chose 
to do the former. The induced steady state polarization 
Pg is defined as the fixed point of the equation ([T|) where 
we assume that the electric field Ef{t) = Eg is constant 
and the polarization P = Pg are time independent. We 
get 

aPg+l3P^^Eg. (5) 

We assume the solution of ([1]) to be of the form 

E = Eg + SE, P^Pg + SP, (6) 

where SE and 5P are small compared respectively to Eg 
Pg. Plugging this into the equations ([T]) we get the lin- 
earized equations around the stationary solution {Eg, Pg) 



n^d^SE d^SE I d^6P 



t2 d^5P -iT^ d6P 
1- _! 

eo df^ eo dt 



eoc" 



dt^ 



5{z) 



(7) 



{a + i(3P'^)5P = 5E{t). 



B. Scattering of linear waves by the thin film 

We now proceed to solve the linearized equations ([7]) 
by using a scattering theory formalism. We separate time 
and space by assuming a periodic solution 



5E = e{z)e' 



5P — pe 



We get 



n^k'^e+ — = k^6{z)p, 



dz 



(8) 



(9) 



{n^ + ijT^u - t'^u?)p = eoe(O), 
where the wave number k — lo/c and where we introduced 
^2 = eo(a-|-3/3p2). (10) 
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In the scattering we assume the electromagnetic wave 
to be incident from the left of the film located at z = 0. 
We then have 



Re' 



,z<0 



e = Te" 



,z > 



(11) 



where R is the amplitude of the reflected wave and T the 
amplitude of the transmitted wave. These expressions 
are related to the A, B, C parameters of the previous 
section through the relations R = B/A, T = CIA. We 
have the following interface conditions at z = [l^, [l^l 



e(0-) = e(0+), 
They imply 

l + i? = T. nT 



Co 



(12) 



(1-i?) 



-ik — p 



Using the second relation of ([9]) to obtain p we get the 
transmission coefficient T 



T = 



2c(»^ + ijT^UJ - T^UJ^) 

c{n + l)(r2^ + ijT^Lj — T^ixP-) + iluj 
The refraction coefficient is 

c(l — n)(f7^ + i^T^ijj — T^w^) — ilijj 



R = 



c(l + n){^^^ + ijt'^uj — r^oj^) + iluj 



(13) 



(14) 



The bound states are the poles of the reflexion and trans- 
mission coefficients. Their existence indicates that the 
system has resonant modes that can be excited by an in- 
coming wave. The real part of the bound states is the 
oscillation frequency and the imaginary part is the in- 
verse of the decay time of the mode. 
The poles of R, T are given by 



c{n + l)(ri^ -I- ijT^Lo - T^w^) + iliu : 
which is the second degree equation 
I 



0, 



2 2 
W T 



{n + l)c 



= 0, 



(15) 



whose roots are 
I 



w = — , , 

2 \{n + l)c 



/f^2 1 

+ 7) 



I 



■1 



[n + l)c 

The imaginary part and real part of w = + i^i give re- 
spectively the decay time Tdec and the oscillation period 
T 



TqsC — 



27r 



(17) 



We will estimate these parameters for the experiment in 
the next section. 



C. Green function solution of the linearized 
equation 



The linearized equation 
t2 d'^SP d5P 



eo dt 



[a -f 3I3P^)6P = 5E{t), 



can be solved to obtain the polarization response 5P{t) to 
a give incoming electric field 5E[t). For that we introduce 
the Green function G{t — to) which satisfies 

r'^+ + ^'G = eo6{t - t,). (18) 

Using the Laplace transform 

G{s) EE / e-'*G{t)dt, 
Jo 

and assuming that G(0) = and dG/dt{Q) = we obtain 



G(,s) 



(19) 



To get the inverse Laplace transform one expands this 
rational function as 



Gis) = 



2(si - S2) \s - si 



-sto 



-sto 



S2 



where Si, S2 are the roots of the denominator of ([T 
This yields the Green function 



eo 



Git - to) - „ , , 
2(si - S2) 

Git - to) = 0, 



gSi(t-to) _ gSi(t-to) 



, t>to 



t < t 







The roots are complex conjugate si = uji + iiOr , §2 ~ 
Ljj — iLUr SO we obtain the final result 



G(t-^o)=e--* ^^"["f t>to 



2uJr 



(20) 



Git - to) - 



t < to 



The polarization response 5Pit) to a given perturbation 
of the electric field 6Eit) is the convolution integral 



6Pit) 



Git - to)5E{to)dto. 



(21) 



D. Induced polarization caused by a short electric 
pulse 

Now let us consider an alternative method to inves- 
tigate the polarization response. Suppose that a ferro- 
electric film is in a polarized state Ps caused by a con- 
stant electric field Eg. At certain moment we send in a 
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short electric pulse so the polarization of the film changes 
during a short time period. After that the polarization 
relaxes to the steady state position Pg. If the electric 
pulse is sufficiently short i.e. close to a " (5 function" then 
the polarization response will follow the Green function 
([20|) . We term this action " 5- function-like pushing" the 
nonequilibrium polarization. 

The linearization of the equation for the polarization 
dl]) near the steady state Ps assuming P ~ Pg + p with 
the condition p « Pg results in the equation for the 
Green function (ITOl) . Let us suppose that the extremely 
short electric pulse acts at t = 0. Then we can conclude 
that the evolution of p{t) after t = is described by 

2uir 

Hence the polarization decay rate is defined by the imag- 
inary part of the complex frequency (|16l) , uji and the cor- 
responding time for the polarization to attain the equi- 
librium is 

Tdec ~ ^I'^i- 

The real part of the frequency (|16l) corresponds to oscil- 
lations of the polarization as it is approaches the steady 
state value. Here we recover the results obtained using 
the scattering formalism. 

III. EXPERIMENTAL RESULTS 

The experiment was performed u sing the nonlinear op- 
tical stroboscopic technique as in [14]. For the second 
harmonic (SH) generation, the radiation of a titanium- 
sapphire laser (MaiTai, New-Port-SpectraPhysics) was 
used with a pulse duration of 100 fs, a wavelength of 
780 nm, a repetition rate of 100 MHz and an aver- 
age power of 100 W. The experiment was performed 
at room temperature. The 70nm thick BaojSrQ ^TiOs 
films were deposited onto a MgO substrate by RF sput- 
tering. For such a composition, the Curie temperature 
equals Tc = 20 C. However for thin films the phase tran- 
sition is blurred around this value and this is confirmed 
by the presence of a narrow hysteresis above Tc [l3]- A 
voltage pulse of duration about 25 ns, produced by an 
Avtech pulse generator, was applied to the copper con- 
tacts on the EST film. The polarization response of the 
ferroelectric film was measured as the coherent SH inten- 
sity [H. 

Figures [H and |3] show the polarization response as a 
function of time of the EST thin film to electric pulses 
having the same amplitudes and opposite polarities. The 
left panels show the incoming electric field pulses as a 
function of time and the right panels show the SH inten- 
sity (proportional to the square of the polarization) as a 
function of time. In figured] the electric field pulse (left 
panel in the figure) is realized by a rapid spiking to a 
zero value from a constant negative voltage background 




time, ns time delay, ns 

FIG. 2: Plots of a 'zero' electric pulse as a function of time 
(left panel) and the subsequent polarization oscillations of the 
film as a function of time (right panel) . The constant field Es 
is negative. 




time,ns time delay, ns 

FIG. 3: Plots of a 'zero' electric pulse and the subsequent 
polarization oscillations of the film as in Figure [2] except that 
the constant field Es is positive. 



with return to the same constant negative value. For fig- 
ure [3] the electric pulse (left panel at the figure) drops 
from its constant positive value to zero and returns to 
the constant. In both measurements the spiking pulse 
will be called the 'zero' pulse. The pulse profile is in- 
tended to have a narrow bell shape, however because of 
certain setup drawbacks a low-amplitude tail appears. In 
the right panels of Figures [5] and [3] it can be seen that 
the polarization oscillates long after the 'zero' electric 
pulse has passed through the thin film. Then the eigen- 
frequency of the polarized thin film can be determined. 

In both pictures the polarization oscillates around its 
stationary values (defined by the constant background 
electric field) with very close oscillations periods, about 
60 ns. This agrees with the relation (1^^ in the limit 
when the second term can be neglected in the square 
root because fl only depends on the amplitude \Eg\. 

The next figure ^ shows the polarization dynamics 
when there is no constant electric field Eg after a delta 
function like electric pulse of negative polarity (shown in 
the left panel) passes through the film. This type of elec- 
tric pulse is the analog of a 'normal pulse' studied exper- 
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20 40 




time, ns 



20 40 

time delay, ns 



FIG. 4: Left panel, plot of the electric field pulse as a function 
of time for a "normal pulse" (see text). The right panel shows 
the subsequent polarization of the thin film as a function of 
time. 



Tc, [C] 


a, [m/F] 
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2.2 • lO'' 


2.52 • 10** 



TABLE L Material parameters for Bao.rSro.sTiOs 



imentally in [lj|. When the film is not polarized the SH 
cannot be generated, so that its intensity is about zero. 
The SH signal from the perturbed unpolarized state is a 
few times smaller than the one for the previous experi- 
ments. There the steady polarized state of the thin film 
was studied by a 'zero' pulse. There are almost no polar- 
ization oscillations when perturbing a non polarized thin 
film. Again this agrees with the estimate 



A. Discussion of the experimental results 

There is a relatively broad range of experimentally ob- 
tained values of the Landau coefficients in EST solid so- 
lutions. The numbers vary with the fabrication method. 
The Landau coefficients found by [l^l for a EST solid 
solution Bao jSrQ^TiO^ in SI units are presented in Ta- 
ble I. In the experiment performed the film thickness is 

1 — 70nm (as in 114]). We assume that the substrate 
refraction index is n = 1.5. These parameters enable to 
calculate the frequency fl. For a normal pulse the ap- 
plied voltage Vs = so Eg = and Pg = 0- This gives 
il^ — ega « 210"^. In the presence of an electric field 
Eg = Vs/l where Vg = lOVolts we have Eg w lO^V/m 
so that from ^ we get Pg « 0.79C/m'^. This gives 

02 ^ gpQ! « 1.810-3. 

First let us estimate the term l[{n + l)c]^^ in the for- 
mula (dD). We have l[{n + l)c\~^ w 10-^*^ which is 
very small so that this radiative damping can be com- 
pletely neglected for this particular experimental situa- 
tion. /^From the experimental data it can be seen that the 
polarisation response for a voltage = is qualitatively 
different from the one for Vg = 10 volts. In particular it 
has a smaller decay time and practically no oscillations. 
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TABLE 11: The decay time Tdec as a function of the parame- 
ters (r, 7) for the zero pulse V3 — volts. 
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TABLE 111: The oscillation period Toac as a function of the 
parameters (r, 7) for the non zero pulse Vs = 10 volts. 



Let us estimate the decay time T^ec from (fTO)) . We have 

Tdec = 2/7 
if (ri^/r^ - 7^/4) > 0, otherwise 
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r 



Tdec = 7: + \l T 



This decay time is shown in Table II for different values of 
the free parameters (r, 7). Clearly if (il^/r^ — 7'^/4) > 
there will be oscillations in the polarization. Comparing 
the estimates of the table with the result shown in Fig. 
|4] indicates that 7 « lO^s-i. This agrees with the value 
obtained in a previous study by the authors [23*1. 

Let us now compare the oscillation period obtained for 
the non zero pulse Vg — 10 volts with the estimate p^ . 
We have 



To. 




(22) 



This value is computed for different values of 7 and t 
and reported in Table III. The zero entries correspond to 
the situation where il"^ /t"^ — j'^ /A < 0. As can be seen 
the values closest to what is observed in the experimental 
plots Figs. [2] and [3] are r = lO'^^s and 7 = 10^. This 
value of the oscillation time r agrees with the one that 
can be computed using the speed of sound Cg — 510'^m/s 
in the material and the film thickness I. We get r = 
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15 30 45 to t5 20 



t(10''s| l(10''s) 

FIG. 5: The polarization response of the film for a normal 
pulse {Es = 0) (left panel) and for a zero pulse (right panel). 

l/cs ~ lO^^^s. The estimate for the bulk material using 
the Newton equation describing the lattice oscillations 
T^j: w Meo/iNe^) « IQ-^^s'"^ where M is the atomic 
mass of BST gives t « 10^ ^'^s which is much too small. 

We have used the expression ([2T|) to compute the po- 
larization response of the system to a given perturbation 
6E of the electric field. We assume a gaussian 6E 

5E{t)=e^p(~^], (23) 

and normalize times by T = 10ns. We chose 7 = 
lO^s^^, T = lO^^'^s. In the normalized units we have 
7 — 0.1, We — 1. The normalized frequency w — 0.45 
for the zero pulse and cj = 4.5 for the normal pulse. We 
calculated the polarization response by numerical inte- 
gration using the trapeze method. To compare with the 
experimental data we computed SP^. The results are 
presented in Fig. [S] for a normal pulse (Eg — 0) in the 
left panel and for a zero pulse in the right panel. One can 
see on the left side of the plots the first response of the 
polarization. We have omitted it because it is the forced 
response due to the probing pulse 5E. We only present 
the subsequent free evolution of SP. For Vg = the po- 
larisation decays with few oscillations. For Vs = lOvolts 
the polarization oscillates much longer. This is in quan- 
titative agreement with the experimental plots presented 
above. 



IV. CONCLUSION 

We analyzed the polarization oscillations occurring in 
a thin ferroelectric film as a short electric pulse crosses 



it. We consider that the ferroelectric is in the para- 
electric phase (high temperature). As in a pump probe 
experiment the film is initially in a static polarization 
state Pg induced by a constant voltage Vg. Using the 
Landau-Ginzburg-Devonshire theory we computed this 
static polarization and the subsequent oscillations of the 
polarization induced by a short voltage pulse. These 
were analyzed using a scattering theory formalism and 
a Green's function approach. Two channels of dissipa- 
tion were identified, a radiative damping and an intrinsic 
damping. 

This theory was applied to explain the experimental 
time evolution of the polarization for a thin ferroelec- 
tric film of BST. The polarisation is estimated indirectly 
through second harmonic generation. For this experi- 
mental situation we show that the radiative damping can 
be neglected and only the intrinsic damping should be 
considered. From a comparison of our model to the ex- 
perimental plots we estimated the important parameters 
T the response time and 7 the relaxation coefficient. Us- 
ing these values our theoretical estimates of the decay 
time and of the oscillation period agree well with the ob- 
servations. In particular for a normal pulse for which 
Vs = the damping dominates and we only see a few os- 
cillations of the polarization. On the contrary for a zero 
pulse for which Vg — lOvolts the damping time is longer 
compared with the oscillation period. Here we obtain 
many oscillations of the polarization. 

Although the radiative damping appears as artificial 
in this particular high temperature situation, at low tem- 
perature phonons become frozen and the radiative damp- 
ing may become predominant. The model can be further 
elaborated by including the spatial inhomogeneity of the 
thin film, the depolarization effects of the boundaries, 
and taking into account the internal stresses of the film. 
Also the dynamics of the polarization of a thin film multi- 
layered structure can be investigated by generalizing the 
model. 
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